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MEIHQ.D 


The  theoretical  basis  of  the  program  is  contained  in  Volume  I 
and  is  not  repeated  in  this  manual. 

Variable  dimensions  of  the  grid  have  been  used  throughout  so 
that  any  size  problem  can  be  run  if  sufficient  computer  size 
and  time  are  available. 

The  solution  scheme  is  a  Gaussian  reduction.  This  gives  good 
answers  with  reasonably  short  run  times. 

The  program  may  be  run  as  a  standard  ax (symmetric  program 
giving  either  one  or  both  of  the  axisymmetric  or  torsion 
solutions.  Multiple  passes  may  be  made  through  the  program 
inputting  the  proper  load  coefficients  to  gain  an  asymmetric 
loading  solution.  The  accumulation  routine  will  accumulate 
the  displacements  and  strains  for  all  harmonics  to  give  the 
final  solution. 

The  program  Is  a  finite-element  program.  To  solve  for  the 
displacements,  stresses  and  strains  throughout  a  body.  It  Is 
divided  into  many  small  pieces  called  elements  by  defining  a 
grid  network  of  nodes  over  a  cross  section  of  the  body.  For 
each  element,  a  stiffness  matrix  Is  defined  according  to  the 
basic  assumptions  for  that  element.  A  force  vector  is  created 
by  applying  the  pressures  and  loads  to  the  nodes.  The  element 
stiffness  matrices  and  force  vectors  are  then  assembled  into  a 
master  set  of  equations,  and  the  equations  are  solved  using  a 
Gaussian  reduction  technique.  This  method  of  solution  has 
been  extensively  used  at  Thlokol  for  many  years  and  has 
yielded  good  answers  in  reasonably  short  run  times. 

There  are  two  basic  elements  available  in  this  program,  and 
one  or  the  other  must  be  selected  in  the  RUN  input  section. 
The  first  is  the  linear  displacement  element.  The  assumption 
for  it  is  that  the  displacements  are  a  linear  function  through 
the  element.  This  element  is  used  for  the  basic  axisymmetric 
runs,  the  asymmetric  loading  runs  and  the  torsion  runs. 

The  second  element  is  the  isoparametric  element.  It  is  a 
quadratic  element  where  the  displacements  are  assumed  to  be  a 
quadratic  function  through  the  element.  This  element  will  In 
general  give  better  results  for  the  same  size  or  a  smaller 
grid.  It  has  not,  however,  been  as  extensively  used  as  the 
1 (near  element. 
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Type  of  Material  section,  there  are  five  element  types 
Types  1,  2,  4  and  5  are  all  linear  elements  and  differ 
the  number  of  degrees  of  freedom  they  can  handle. 


LQM  DEFINITIONS 


For  the  single-pass,  axlsymmetrlc  run,  loads  are  Input 
directly  as  they  apply  to  the  entire  body  of  revolution. 

For  asymmetric  loading,  a  Fourier  expansion  of  the  form 


N  M 

f(fl)  ■  A0  ♦  Z  Afl  cos  n0  +  Z  B  sin  n 0 

n*l  n»l 

must  be  generated,  where  N  is  the  number  of  terms  needed  to 
accurately  represent  the  load.  All  loads  are  defined  In  terms 
of  this  function  of  theta.  The  loads  are  then  Input  with  the 
constant  term  Ap  In  the  first  pass  and  either  Ap  or  Bp  being 
entered  for  the  n+lst  pass. 

The  assumption  is  Implicit  In  the  program  that  the  loads  are 
symmetric  about  the  0°  -  IS 0°  plane  through  the  body.  Because 
of  this  symmetry,  the  radial  and  axial  loads  can  go  In  only  as 
an  even  function  and  the  tangential  loads  only  as  an  odd 
function.  That  Is,  only  the  A  terms  can  be  used  for  R  or  Z 
loads,  and  only  the  B  terms  can  be  used  for  the  e  loads. 

The  units  of  the  function  f  and  the  coefficients  A  and  B  must 
bo  in  pounds  per  square  inch  for  pressure,  or  pounds  for  nodal 
point  load,  over  the  entire  surface  affected. 

Pressure,  shear  and  traction  loads  act  on  the  entire  surface 
of  revolution  on  which  they  are  applied.  That  Is,  the  total 
load  is  equal  to  x DLP  where  D  Is  the  diameter,  L  the  length, 
and  P  the  pressure  or  shear. 

Nodal  point  forces  apply  a  load  directly  to  a  point  and  are 
not  distributed  over  any  surface.  They  are  a  load  on  a  ring 
as  the  loaded  point  is  revolved  about  the  axis. 

See  Appendix  B  for  the  description  of  a  program  to  generate 
the  Fourier  coefficients. 
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PROGRAM  COORDINATES 


The  normal  orientation  of  the  coordinate  system  Is  with  the 
R-axis  horizontal  and  the  Z-axis  vertical.  This  is  referred 
to  here  as  a  right-handed  system;  that  Is,  R  rotated  Into  Z 
gives  a  vector  out  of  the  paper. 

The  indexing  system  of  the  program  can  also  be  layed  out  as  a 
coordinate  system.  By  definition,  a  right-handed  system  is 
one  in  which  the  I  axis  rotated  Into  the  J  axis  gives  a  vector 
out  of  the  paper. 

The  two  systems,  coordinate  and  Indexing,  must  always  agree. 
That  is,  they  must  both  be  either  right-handed  or  left-handed. 
If  one  is  right  and  the  other  left,  the  program  will  give 
negative  area  errors  on  the  grid  and  stop  execution. 

All  angles  in  the  system  are  measured  from  the  positive  R  axis 
counterclockwise  in  the  right-handed  system.  If  the  system  Is 
transformed,  this  relative  orientation  must  be  maintained. 

There  are  "local  origins"  referred  to  in  the  input  section. 
These  are  simply  temporary  origins  to  aid  in  the  designating 
of  that  one  node  or  line.  The  point  will  be  translated  back 
into  the  global  coordinate  system  by  adding  the  coordinates  of 
the  local  origin  to  the  coordinates  of  the  point. 


lilt  illit  GENERATOR 


By  means  of  the  Line  Generator  Record,  line  and  arc  segments 
may  be  generated  Internally  to  connect  any*  two  points  with  a 
set  of  nodes.  The  following  options  are  available: 

A.  Equal  interval  -  In  the  straight  line  all  intervals 
between  nodes  are  equal.  In  an  arc  the  angular 
intervals  are  equal. 

B.  Square  root  of  r  -  Where  M  Is  the  total  number  of 
Intervals  and  I  goes  from  1  to  fl-1,  each  Interval  Is 
(2i-l)/N2  of  the  total  length.  For  a  line,  this  Is 
an  increment  of  length;  for  an  arc,  it  I s  an  angular 
Interval . 

C.  Geometric  progression  -  The  increments  of  length  or 
angle  will  be  In  a  geometric  progression.  Any  of  r 
(the  common  ratio),  a  (the  first  Interval),  or  l  (the 
last  interval)  may  be  specified.  In  addition,  the 
second  node  in  the  line  may  be  specified  and  r  will 
be  calculated  using  the  interval  from  the  first  node 
to  the  second  node  of  the  line  as  a. 

Any  of  these  options  may  be  used  to  generate  any  line 
in  the  figure.  The  lines  may  go  In  either  ascending 
or  descending  order  of  I  or  J. 


♦Either  I  or  J  must  be  constant  for  all  the  nodes  In  a  line  or 
arc. 
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Iii£  fifilU  generator 


The  grid  generator  used  In  this  program  Is  one  developed  at 
Thlokol  by  Dr.  William  Cook.*  This  is  a  linear  Interpolation 
method  which  yields  a  good  grid  under  most  conditions 
regardless  of  the  curvature  of  the  boundaries. 

A  two-dimensional  space  is  defined  in  5  and  ij ,  and  the  grid 
boundary  nodes  are  napped  Into  this  space.  A  linear 
interpolation  scheme  Is  then  used  to  fill  in  the  remainder  of 
the  £  ,  tj  space.  Two  functions,  f  and  g,  are  then  defined 
that  map  t lie  C  ,  R  space  Into  the  X,Y  space  of  the  grid.  These 
functions  are  then  used  to  define  the  remainder  of  the  grid. 

There  may  be  any  number  of  major  partitions  in  the  grid  and 
any  number  of  subpartitions  within  each  of  these. 

A  major  partition  Is  defined  as  a  closed  rectangular  section 
of  matrix  (the  figure  may  be  of  any  shape,  but  the  section  of 
matrix  must  be  capable  of  being  defined  by  a  minimum  and 
maximum  of  I  and  J),  ail  of  which  will  be  generated  by  one  and 
only  one  grid  generation  record.  The  subpartitions  are  closed 
rectangular  sections  of  matrix  lying  within  the  major 
partition. 

All  boundaries  of  each  partition  must  be  completely  specified, 
and  there  must  be  no  nodes  specified  that  do  not  lie  on  the 
boundaries  of  a  partition. 

Each  major  partition  must  be  specified  on  a  Grid  Generator 
Record,  but  the  program  will  pick  up  all  subpartitions  lying 
within  the  specified  major  partition. 

Due  to  the  nature  of  the  Cook  grid  generator,  the  following 
restrictions  must  be  observed  when  using  it: 

1.  All  sections  of  the  matrix  that  are  to  be  generated 
by  it  must  be  rectangular  and  the  boundaries  must  be 
completely  defined. 

2.  If  interior  nodes  are  specified,  they  must  form 
rectangular  subsections  of  the  section  of  matrix 
being  generated.  There  is  no  limit  to  the  number  of 
such  subsections  allowed. 


*Cook,  William  A.,  BODY  ORIENTED  (NATURAL)  CO-ORDIMATES  FOR 
GENERATING  THREE-DIMENSIONAL  MESHES,  International  Journal  for 
Numerical  Methods  in  Engineering,  Vol  3,  1974,  pp.  27-43. 
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INPUTTING  liiL  PROGRAM 


There  are  several  things  about  the  program  that  you  should  be 
aware  of  as  you  set  up  your  Input.  The  sections  below  define 
both  general  things  to  keep  in  mind  and  the  specific  methods 
of  input  required  for  the  various  options.  For  additional 
information,  see  Vol  It  of  this  report. 

i  General 

The  RUN  section  must  always  be  used  to  specify  both 
the  geometry  type  and  the  element  type.  In  this 
version  of  the  program,  the  geometry  type  is  always 
axi symmetric  (option  3).  The  element  type  can  be 
either  linear  displacement  (option  4)  or 
isoparametric  (option  5). 

Each  node  defines  an  element,  except  for 
isoparumeterlc,  where  it  is  every  other  node.  The 
elements  are  referenced  by  the  node  with  the  smallest 
indices  of  the  nodes  comprising  the  element.  Those 
nodes  which  lie  on  the  boundary  of  the  geometry  and 
do  not  reference  a  real  element  are  given  a  type  code 
of  9,  which  signifies  a  missing  element.  This  type 
code  of  9  is  also  used  to  code  an  internal  element 
which  has  no  material  in  it.  Type  9  elements  require 
no  material  code,  and  this  can  be  set  to  0  or  left 
blank. 

Plots  can  be  produced  quite  simply  by  specifying  the 
orientation  desired,  the  range  of  indices  of  the 
nodes  to  be  plotted,  and  the  Y  paper  size.  The  data 
will  then  be  scaled  to  fit  the  paper  height,  and  the 
same  scale  factor  will  be  used  for  the  X  direction. 
The  actual  scale  values  used,  the  minimum  values  and 
the  maximum  values,  if  calculated,  will  be  printed 
when  each  plot  is  produced. 

Several  Input  sections  use  Index  increments  referred 
to  as  II  and  Jl .  The  value  II  is  an  I  index 
increment,  that  Is,  we  go  from  II  to  12  In  increments 
of  II.  The  value  Jl  is  a  J  index  increment,  and  we 
go  from  Jl  to  J2  in  Increments  of  Jl.  For  example, 
if  we  have  a  grid  where  we  have  material  1  In  the 
elements  with  an  even  J  Index  and  material  2  In  the 
elements  with  an  odd  J  index,  we  can  apply  these  in 
only  2  records  by  using  a  Ji  of  2.  Another  example 
of  its  use  is  in  the  Isoparametric  input.  This  Is 
explained  more  fully  in  section  VII  below. 


11 


I 


Ih£  Pflsic  Axisvmmetr i c  Problem 


This  Is  the  case  for  which  there  are  no  Iterations 
and  the  basic  linear  displacement  element  Is  being 
used.  The  RUN  section  must  show  this  by  Inputting 
options  3  and  4. 

1 1  Asmafit-rlc  loading 

The  asymmetric  loading  of  an  axl symmetric  body  is 
accomplished  by  making  multiple  passes  through  the 
program  within  one  run.  Each  pass  has  the  harmonic 
coefficients  for  one  term  in  the  Fourier  expansion 
entered  as  load  conditions.  Run  option  9,  asymmetric 
loading,  must  be  specified  in  each  pass. 

The  node  number/  Lll  in  the  General  Data,  and  the 
highest  harmonic  allowed,  L12  in  the  General  Data, 
must  both  be  specified  in  each  pass,  and  the  highest 
harmonic  allowed  must  be  nonzero. 

The  accumulation  of  the  output  is  necessary  for 
proper  interpretation  of  the  results.  Each  pass  will 
give  results  that  are  the  maximum  for  each  coordinate 
direction.  These  values  have  meaning  only  when 
multiplied  by  the  proper  sine  and  cosine  terpis  and 
added  to  the  other  terms  for  that  node  or  element. 

The  undeformed  geometry  can  be  plotted  at  any  time  in 
the  run.  Accumulation  and  plotting  of  the  deformed 
geometry  must  be  in  the  last  pass  only.  Plots  can  be 
made  of  the  displaced  geometry  at  any  angle  0  or 
axial  location  Z  by  accumulating  the  data  where 
desired  and  then  inputting  a  plot  section  immediately 
after.  The  plot  routines  will  always  plot  the  last 
data  accumulated.  See  the  sample  input  for  on 
example  of  this. 

V  Incremental  LPildJaE 

When  tills  option  is  used,  the  program  automatically 
cycles  through  as  many  tines  as  there  are  load 
Increments.  The  cumulative  totals  of  the 
displacements  and  strains  may  be  printed  out  for  each 
increment  under  the  control  of  the  Selective  print 
option.  The  print  suppress  option  (output  option  1) 
and  the  Selective  print  section  control  the  output  of 
the  final  results,  and  the  incremental  loading  print 
flag  with  the  Selective  print  option  controls  all 
others . 
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V 


VI 


All  loads  and  boundary  conditions  may  be 
Increments.  The  fraction  of  the  load  to 
in  each  step  may  be  specified,  or  on 
Increments  will  be  calculated  in  an 
progression. 


appl led  in 
be  appl l ed 
a  flag  the 
ar i thmet 1 c 


A  set  of  flags  is  provided  so  that  specified  loads 
may  not  be  applied  incrementally. 


Each  iteration  is  handled  essentially  as  a  separate 
run  with  the  program  looping  through  the  entire 
solution  scheme.  The  displaced  geometry  reflecting 
all  the  previous  loading  is  used  in  each  Iteration  so 
that  stiffness  due  to  bending  will  be  taken  into 
consideration. 


Large  Deformation 

The  use  of  this  option  takes  into  consideration  the 
nonlinear  terms  In  the  strain  displacement 
relationship  and  iterates  for  the  final  solution. 
You  can  also  enter  strain  dependent  material 
properties  to  take  into  consideration  the  nonlinear 
material  behavior. 


A  Convergence  criteria  and  a  series  of 
underrelaxation  factors  must  be  entered.  The 
underrelaxation  factors  must  be  positive  values  no 
greater  than  1.0,  and  at  least  the  last  two  values 
must  be  1.0.  The  program  will  iterate  until 
convergence  or  the  maximum  number  of  iterations  is 
sat i sf ied. 

As  In  the  incremental  loading  option,  the  print 
suppress  flag  controls  the  final  output  and  the  large 
deformation  print  flag  controls  nil  others. 

The  energy  calculation,  reaction  forces  and  accuracy 
check  are  not  valid  for  this  option  and  ore  not 
a  1 1  owed . 


Incremental  Load  I  or.  SM  Large  Deformation 

These  two  options  may  be  combined  to  allow  a  large 
deformation  Iteration  inside  each  incremental  loading 
loop.  The  controls  on  each  remain  the  same,  and  the 
same  set  of  underrel axat i on  factors  will  be  used  for 
each  increment  of  load. 
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VII  Isoparametric  Elements 

The  i soparametr ic  element  in  this  program  is  a 
quadratic  element  with  three  nodes  on  each  side.  The 
geometry  is  input  and  generated  just  ns  if  a  linear 
displacement  element  were  being  used.  Note,  however, 
that  all  element  boundaries  and  naterial  boundaries 
must  fall  on  nodes  with  odd  Indices. 

The  internal  handling  of  this  element  requires  some 
special  consideration  in  the  input.  The  element  type 
codes  and  material  numbers  must  only  be  on  the  nodes 
which  designate  element  corners.  To  achieve  this,  II 
and  Jl  in  the  TYPE  section  must  be  input  as  2's. 

VI  I  I  Torsion 

The  torsion  solution  is  part  of  the  zero  harmonic  in 
the  asymmetric  loading  case,  because  it  is  the  only 
element  at  present  that  has  a  theta  degree  of 
freedom.  A  RUN  parameter  of  9  must  then  be 
specified. 

There  are  to  different  ways  that  torsion  can  be  run: 

(1)  with  the  full  asymmetric  calculation  givln<;  all 
three  degrees  of  freedom,  and  using  a  type  4  element: 

(2)  with  torsion  only,  in  which  case  a  type  5  element 
will  be  used.  The  type  5  element  Is  faster  than  the 
type  4.  Accumulation  can  be  made  for  both  elements, 
or  for  any  combination  of  2,  4  or  5  elements,  since 
they  are  all  asymmetric.  Note  that  the  torsion  loads 
act  on  a  moment  arm  equal  to  the  radial  coordinate  of 
the  point  where  they  are  applied. 

IX  StflbM  i  tY 

This  option  is  used  to  determine  the  stability  of  a 
flat,  vertically  stacked  bearing.  The  only  input 
sections  required  are  the  General  Input,  the  Run 
input,  and  the  Stability  input.  The  values  in  the 
General  section  will  not  be  used,  but  due  to  program 
flow  the  section  must  be  included.  Run  option  20 
must  be  specified. 

X  Service  Li fe 

To  determine  the  service  life  of  a  bearing,  this 
program  must  be  used  in  conjunction  with  program 
S3359SL.  When  a  run  is  made  where  data  needs  to  be 
saved  for  the  service  life  calculations,  output 
option  20  must  be  specified.  The  stresses  and 
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strains  will  then  be  output  on  a  fiie  v/hose  name  is 
SERVLIFE.  Each  set  of  data  put  on  this  file  will  be 
identified  by  two  Integers  which  are  printed  on  the 
output:  the  first  Is  the  julian  date  and  the  second 
is  the  time  of  day  in  cent l seconds.  These  will  be 
used  to  identify  the  set  of  data  for  program  S3359SL. 
These  sets  can  be  put  on  separate  reels  of  tape,  or 
they  ran  be  stacked  together  on  one  tape  by  using  the 
proper  Job  Control  Language. 


TIMING 


Timing  on  this  program  Is  a  problem  because  so  many  factors 
affect  the  execution  speed. 

In  general  the  time  required  for  a  run  is  a  function  of  I2  and 
J. 

The  asymmetric  loading  problem  requires  about  two  times  as 
long  for  each  pass  as  a  simple  axisymmetric  problem. 

Each  iteration  on  iterative  runs  must  be  counted  as  a  full 
execution. 

From  our  experience,  the  best  way  to  estimate  the  time  is  by 
experience  with  the  guidelines  given  above.  However,  the 
following  equation  will  give  a  starting  point: 

Time  (in  minutes)  ■  F  I2  J  +  1.2 


where  F  ■  .0012  for  an  asymmetric  run  and  .0006  for  an 
axisymmetric  run,  I  is  the  I -dimens ion  of  the  grid,  and  J  is 
the  J-dimension  of  the  grid.  This  equation  is  for  the  IBM 
370/155;  the  factor  F  and  the  constant  1.2  may  change  for 
other  machines. 
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RESTRICTIONS 


1. 

The  isotropic  Poisson's 

ratio 

is 

restricted  to 

0  <  v 

<  .5.  At 

least  under 

some 

ci rcumstances 

a  zero 

Pol sson 's 
s 1 ngular . 

ratio  will 

make 

the 

stiffness 

matrix 

2. 

1 f  a  node 

is  on  the 

axi  s 

of 

revolution 

of  an 

ax  I  symmetrical ly  loaded  body,  the  radial  displacement 
of  that  node  must  be  fixed  at  zero. 


3.  At  least  one  axial  component  of  displacement  must  be 
specified  to  avoid  rigid-body  axial  motion.  A 
sliding  boundary  condition  may  also  fulfill  this 
requirement  if  It  restrains  axial  motion. 

4.  In  any  row  of  the  grid,  I  must  run  consecutively 
through  each  value  IMIN  £  I  i  I  MAX  for  that  row. 
Type  9  elements  may  be  used  where  needed  to  provide 
“voids"  In  the  grid. 

5.  In  the  input,  a  maximum  of  9  digits  is  allowed  on 
each  side  of  the  decimal  point  in  a  number. 

6.  There  is  a  maximum  of  32,767  materials  allowed  In  the 

program.  However,  due  to  the  special  internal 

handling,  the  number  of  materials  used  with 

isoparametric  elements  is  limited  to  50. 
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VARIABLE  DIMENSIONS 


To  facilitate  running  of  both  large  and  small  problems#  the 
working  arrays  In  the  program  have  been  variably  dimensioned. 
Two  things  should  be  noted  here  about  the  program.  First#  the 
program  will  take  as  much  core  storage  as  It  is  given;  that 
Is#  on  the  step  resource  usage  given  at  the  end  of  the  run# 
the  amount  of  core  used  and  requested  wi 1 1  always  be  the  same. 
Second#  the  actual  amount  of  working  storage  used  in  each 
phase  will  be  printed  out  as  part  of  the  output. 

The  table  on  the  following  pages  gives  the  maximum  core 
required  for  a  run.  There  are  many  factors  which  affect  the 
core  required  for  any  given  run#  so  the  first  guess  should  be 
taken  from  the  table  and  then  adjusted  after  a  complete  run 
has  been  made.  The  excess  core  available  for  each  part  of  the 
program  is  given  with  the  output. 

To  find  the  region  required  for  a  run#  a  look  in  the  table  for 
the  I  and  J  size  of  the  grid  will  give  the  region  in  K  bytes. 
If  an  asymmetric  loading  case#  or  an  isoparametric  case  is  not 
being  run#  the  region  may  be  reduced  by  up  to  one-third. 

For  a  stability  calculation#  only  190K  will  be  needed  for  the 
run. 
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?as 


Emm  input 


Input  to  the  program  Is  via  subroutine  FREFRM.  This  is  a 
record-oriented  free-form  Input  routine.  It  allows  data  to  be 
entered  without  regard  to  card  columns  and  yet  be  input  as 
individual  unique  records. 


This  routine  has  the  following  characteristics: 


1.  Numbers  are  not  restricted  to  any  particular  columns 
on  the  card. 

2.  The  numbers  may  be  separated  with  any  of  the 
fol lowing: 

a.  A  comma  is  a  separator. 

b.  The  sign  of  the  number  is  a  separator. 

c.  The  single  quotation  mark  (')  when  defining  an 
alpha  string  is  also  a  separator. 

d.  The  end  of  a  card  (column  72)  is  a  separator; 

that  is,  if  no  other  separator  Is  encountered 
before  the  end  of  a  card,  that  end  will 

terminate  the  number.  Note  that  an  end  of  card 
in  an  alpha  string  is  an  error. 

e.  The  L  of  an  L-number  is  a  separator. 

f.  A  record  terminator  (;)  will  terminate  the 

number  as  well  as  the  record.  A  record 

terminator  in  an  alpha  string  is  an  error.  A 
string  of  dissimilar  separators  is  only  one 
separator.  For  example,  a  comma  followed  by  an 
end  of  card,  followed  by  an  L  or  quote.  Is  still 
only  one  separator. 

3.  L-Numbers  may  be  used  to  specify  the  relative 
location  of  the  value  in  the  record.  An  L-Number  is 
an  integer  preceded  by  an  L  which  sets  the  counter  to 
a  specific  location  In  the  record,  the  first  location 
being  1.  Thus  the  third  vocation  is  given  by  L3. 
The  L-Number  must  be  separated  from  the  following 
value  by  a  valid  separator.  The  number  must  always 
follow  the  L  immediately.  The  presence  of  anything 
other  than  a  digit  immediately  following  the  L  will 
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cause  an  error 


A  logical  record  must  be  terminated  by  a  semicolon 

(;). 


Alphanumeric  data  can  be  entered  anywhere  in  the 
string  but  must  be  set  off  with  single  quotation 
marks  (').  The  quotation  mark  is  a  separator.  Two 
consecutive  quotation  marks  will  cause  a  quotation 
mark  to  be  entered  in  the  string.  There  are  8 
characters  in  each  location  as  defined  by  the 
L-Numbers. 


If  two  commas  appear  consecutivel y,  no  value  is  input 
between  them.  The  value  previously  put  in  that 
location  will  remain  there. 


If  two  signs  appear  consecutively,  a  zero  is  assumed 
between  them. 


If  a  record  begins  with  an  alpha  character.  It  will 
be  assumed  to  be  a  flag  record  or  a  title,  and  no 
other  data  may  be  entered  in  it.  These  records  nay 
begin  with  'L1,  which  usually  signifies  an  L-Number, 
if  and  only  if  the  second  character  in  the  string  is 
al phabet ic. 


Comments  may  be  punched  in  any  card  by  punching  an 
asterisk  (*)  followed  by  any  desired  comments.  The 
asterisk  and  all  data  following  it  in  that  card  will 
be  ignored. 


Only  the  first  72  columns  of  each  card  will  be  used, 
so  columns  73-80  nay  be  used  for  sequence  numbers  if 
des I  red. 


Data  entered  In  a  location  in  a  record  stays  in  that 
location  In  subsequent  records  until  it  is  changed  by 
entering  a  new  value  in  that  location. 

A  maximum  of  J  digits  is  allowed  on  each  side  of  the 
decimal  point  In  a  number. 


FLAG  AND  TITLE  RECORDS 


Any  record  which  begins  with  alphabetic  data  is  considered  to 
be  a  flag  or  title  record.  Flag  records  are  used  to  flag  the 
beginning  of  each  group  of  records.  If  the  first  four  columns 
of  the  record  contain  a  flag  sequence,  the  record  Is 
considered  to  be  a  flag  record.  If  not.  It  Is  considered  to 
be  a  title  record  and  its  use  Is  determined  from  its  location 
in  the  input  stream.  The  record  terminator  must  not  appear  in 
the  string.  A  title  record  may  span  several  cards,  but  each 
flag  or  title  record  must  end  with  a  record  terminator.  A 
flag  record  cannot  exceed  a  single  card. 

In  the  following  pages,  only  four  characters  are  given  for 
each  flag  since  this  is  all  that  is  checked.  Additional 
characters  may  be  added  to  make  a  full  statement  if  desired. 
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INPUT  ££R  PROGRAM 


As  many  of  the  following  groups  as  needed  may  be  entered. 
Some  of  the  groups  are  required;  others  are  optional. 

The  Title  records  and  General  Data  must  be  the  first  two 
groups  input.  All  other  groups  may  be  Input  In  any  order, 
with  the  only  restriction  being  that  of  logical  sequence  of 
events.  Obviously  you  can  not  use  an  Item  of  data  until  It  is 
input;  therefore.  Line  Generators  must  follow  the  Mode  records 
that  define  their  end  points,  etc.  In  general,  any  group  may 
be  input  more  than  once;  however,  this  practice  should  be 
restricted  to  avoid  overly  complex  input  sets. 

Each  case  run  through  the  program  Is  considered  to  be  a 
separate  case,  and  all  Information  must  be  entered  for  each 
case.  This  is  also  true  of  the  various  parts  of  an  asymmetric 
loading  run.  Except  for  the  accumulation,  each  case  Is 
considered  to  be  completely  separate. 

Note  that  only  columns  1  through  72  may  be  used  for  data. 

LLU.S 

As  many  cards  as  desired  may  be  entered  here,  with  the 
last  card  containing  a  record  terminator.  Each  title 
card  must  start  with  an  alphabetic  character. 

General  Data  -  The  flag  Is  GENE 


LI 

Minimum  1  value 

L2 

Minimum  J  value 

L3 

Maximum  1  value 

L4 

Maximum  J  value 

L5-L9 

Not  used  at  present 

L10 

Base  temperature 

Lll 

Mode  number  If  asymmetric  loading 

L12 

Highest  mode  allowed  if  asymmetric  loading 
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Run  Parameters  -  The  Flag  Is  RUN 


Run  options  will  be  selected  from  the  following  table 
by  entering  the  option  numbers  as  a  sequence  of 
integers.  The  blank  options  represent  options  not 
available  in  this  version  of  the  program. 


1. 


2. 

3. 

4. 

5. 

6. 


Axi symmetric  geometry 
Linear  displacement  element 
Isoparameter ic  element 


7. 

8. 


9.  Asymmetric  loading  of  an  axi symmetric  geometry 

10.  flesh  only  -  check  input  and  test  the  grid. 
Geometry  plots  v/ill  be  produced  if  requested. 

11. 


12. 

13. 


14. 

15. 

16. 

17. 

18. 

19. 

20.  Stability  calculation 
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Output  Options  -  The  Flag  is  OUTP 

The  options  are  selected  by  entering  the  option  numbers 
from  the  following  table.  These  options  apply  to  the 
entire  output  set;  if  output  is  desired  for  specified 
sets  of  nodes  or  elements  only,  then  use  the  Selective 
Print  option. 

If  no  output  options  are  specified,  the  output  will  be 
the  displacements  for  each  node  and  the  stresses  and 
strains  for  each  element. 

1.  Print  Suppress.  The  displacements  and  stresses 
will  not  print  except  as  specified  by  the 
Selective  Print  option. 

2.  Punch  the  deformed  grid. 

3.  Print  the  deformed  grid. 

4.  i 

i 

5. 

6. 

7. 

8.  Accuracy  check.  Calculate  (F-KU)  for  each 
element  and  print  the  10  elements  with  the 
highest  error. 

9.  Reaction  forces.  Print  the  reaction  forces  for 
the  nodes  with  displacement  boundary  conditions. 

10. 

11. 

12. 

13. 

14.  Print  the  strain  energy. 

15.  Print  the  element  material  properties. 

16.  Print  the  loads  for  each  element  with  nonzero 
loads  and  the  boundary  conditions  for  each  node 
with  nonzero  boundary  conditions. 
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17. 

18. 

19. 

20.  Output  data  on  tape  for  service  life 
calculation. 


Special  Points  -  The  Flag  Is  SPEC 

This  section  allows  for  a  maximum  of  100  special  points 
that  are  not  part  of  the  grid.  Their  purpose  is  to 
provide  convenient  points  to  be  used  as  local  origins 
for  nodes,  arcs,  etc. 

LI  K  The  special  point  number.  It  will  be 

referenced  by  this  number  whenever  It  is  used. 

L2  X  or  R  coordinate 

L3  Y,e  or  Z  coordinate 

L4  If  nonzero  the  values  in  L.2  and  L3  are 

considered  to  be  R  and  9,  In  degrees.  If  zero, 
they  are  considered  to  be  X  and  Y. 

L5  The  number  of  a  previously  input  special  point 
may  be  Input  here,  and  it  will  be  used  as  a 
local  origin  for  point  K. 


Nodes  -  The  flag  is  NODE 

Each  of  these  records  defines  one  geometry  point  in  the 
grid. 

LI  I 

L2  J 

L3  R ( I , J ) 

L4  Z( I , J)  or  9  in  degrees 

L5  Spherical  flag  -  Input  a  +1  If  R  and  9  are 

input  instead  of  R  and  Z. 
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L6  Local  origin  flag 

0  -  no  local  origin 

1  -  special  point  whose  number  is  in  L7  is 
the  local  origin 

2  -  the  node  whose  indices  are  given  in  L7 
and  L8  will  be  used  as  the  local  origin 

3  -  the  coordinates  of  the  local  origin  are 
given  in  L7  and  L8 

L7  K,  I  L  or  RL  depending  on  L6 

L8  JL  or  ZL  depending  on  LC 

Line  Generators  -  The  flag  is  LIME 

Each  record  generates  a  sequence  of  node  points  between 
the  existing  end  points.  Either  II  *  12  or  J1  *  J2 
must  hold. 


LI 

11 

Beginning  node 

of 

the  line. 

L2 

J1 

L3 

12 

Ending  node  of 

the 

line. 

L4 

J2 

L5  Rotation  code 

+1  if  rotating  from  the  positive  R  to  the 
pos i t  i  ve  Z  axis. 

-1  if  rotating  from  the  positive  Z  to  the 
posi tive  R  axis. 

L6  Nodal  spacing  option  code 

'O'  equal  interval  option 

*1'  the  square  root  of  r  option 

'2'  geometric  progression  with  'r' 
specified 
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'3'  geometric 

specified 

progression 

wi  th 

•a' 

'4'  geometric 
specified 

progression 

wi  th 

'l' 

•  5 '  geometric  progression  with  the  second 
node  In  the  line  defined  previously  and 
used  to  define  'a' 

L7  This  field  will  contain  a,  r  or  l 

depending  on  which  option  was  specified  In 
L6. 

L8  GMIN  -The  smallest  Interval  allowed  for  a 

geometric  progression. 

L9  GMAX  -  The  largest  Interval  allowed  for  a 

geometric  progression. 

L10  If  the  sequence  of  nodes  is  to  form  an 

arc,  the  way  in  which  the  arc  center  is 

given  is  specified  here. 

0  -  not  arc  -  generate  a  straight  line 

1  -  the  special  point  whose  number  is 
given  In  Lll  is  the  arc  center 

2  -  the  node  whose  indices  are  given  In 
Lll  and  L12  is  the  arc  center. 

3  -  the  center  of  the  arc  is  at  the 
coordinates  given  In  Lll  and  L12 

Lll  K,  1C  or  RC  depending  on  L10 

L12  JC  or  ZC  depending  on  L10 

Bsarinfi  Generator  -  The  flag  Is  SEAR 

This  section  contains  parameters  necessary  to 
generate  a  spherical  or  conical  bearing  model  only. 

The  node  points  will  be  generated  and  element  type 
and  material  codes  will  be  set  as  specified  for  the 
elements  generated.  Type  of  Material  records  do  not 
need  to  be  Input  for  the  elements  of  the  bearing. 
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See  Figure  1  for  the  relative  location  of  the 
indices . 


LI 

L2 

L3 

L4-6 

L7 


L8 

L9 

L10 

111 

L12 

L13 

L14 

L15-18 


L15 


The  minimum  I  value  in  the  bearing.  This 
may  be  any  value.  Default  of  1.  IMNB 

The  maximum  I  value  in  the  bearing  IMXB 

The  minimum  J  value  in  the  bearing. 
Default  of  1,  JMMB 

Not  used  at  present. 

I  EQU  If  this  is  nonzero  the  nodes  will  be 
equally  spaced  across  the  bearing  in  the  I 
direction  Instead  of  having  the  first  and 
last  interval  half  the  others. 

Element  type  -  see  Type  of  Material 
section  for  valid  types. 

♦0  if  a  spherical  bearing. 

+1  if  a  conical  bearing. 

The  inside  radius  of  the  bearing  for 
spherical  bearings  only. 

Radial  coordinate  for  node  IMNB/  JIUIB. 

Radial  coordinate  for  node  IMIIB,  JMXB*. 

Radial  coordinate  for  node  IMXB,  JMtlB. 

Radial  coordinate  for  node  IMXB,  JMXB*. 

Contain  the  axial  coordinates  of  the 
corner  nodes.  For  a  spherical  bearing 
their  only  function  is  to  place  the  node 
in  the  proper  quadrant;  they  need  not  be 
accurate  since  they  will  be  recalculated 
to  be  sure  that  the  node  lies  on  the 
proper  circle.  For  a  conical  bearing  they 
wi 1 1  be  used  as  input . 

Axial  coordinate  for  node  IMNB,  JIUJB. 


♦  JMXB  is  Ji i AX  for  the  bearing  and  is  internally  calculated. 

xn 


L16 

Axial 

coordinate 

for 

node 

IMNB, 

JMXB* 

LI  7 

Axial 

coordinate 

for 

node 

IMXB, 

JMNB . 

L18 

Axial 

coordinate 

for 

node 

IMXB, 

JMXB* 

The  complete  geometry  goes  from  IMIN,JMIN  to  IMAX,JMAX.  Parts 
A  and  C  are  attachment  appliances  on  the  ends  of  the  bearing 
B.  The  bearing  Is  defined  between  Indices  IMNB,  JMNB  and 
IMXB,  JMXB. 

L19  TS  thickness  of  the  shim. 


L20  NS  number  of  columns  in  one  shim. 

L21  ISN  number  of  shims. 

L22  TE  thickness  of  the  elastomer. 


L23  NE  number  of  columns  in  one  elastomer. 

L24  IEN  number  of  elastomers. 


L25-29  Not  used  at  present. 


*JMXB  is  JMAX  for  the  bearing  and  is  internally  calculated. 
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L30 

MS  material  number  for  shim. 

Default  of  1 

L31 

ME  material  number  for  elastomer. 

Default  of  2 

L32-49 

Mot 

used  at  present. 

L50-99 

♦♦SHIMT  shim  thickness  for  each  of  the 
shims  If  not  constant.  Any  value  left 

blank  or  zero  will  be  replaced  by  the 
constant  value  In  L19. 

L100-149 

♦♦ELAST  elastomer  thickness  for  each 
elastomer  If  not  constant.  Any  value  left 
blank  or  zero  will  be  replaced  by  the 
constant  value  In  L22. 

L150-199 

♦♦MATS  material  number  for  the  shims  If 
not  constant.  Any  value  left  blank  or 

zero  will  be  replaced  by  the  constant 
value  In  L30. 

L200-249 

♦♦MATE  material  number  for  the  elastomers 
If  not  constant.  Any  value  left  blank  or 
zero  will  be  replaced  by  the  constant 
value  In  L31. 

Grid  Generator 

-  The  flag  Is  GRID 

Ll 

L2 

11 

J1 

These  specify  the  region  to  be 
generated.  The  generator  used  Is 
the  Cook  grid  generator,  and  the 
region  specified  must  be  a  closed 
rectangular  system  In  the  I,  J 
plane. 


L3 

12 

L4 

J2 

♦♦Starting  with  JMNB  and  proceeding  In  order  out  to  JMXB 
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Type  of  Material  -  The  flag  Is  TYPE 


Each  of  these  records  assigns  a  material  number  and 
an  element  type  to  a  block  of  elements.  A  material 
property  table  must  be  input  for  each  number  assigned 
here.  Note  that  these  are  element  indices. 

LI  II 

L2  J1 

L3  12 

L4  J2  j 

L5  Element  type 

The  valid  element  types  are: 

1  Axl symmetric  with  R  and  Z  degrees  of  freedom 

2  Axl symmetric  with  R,  9  and  Z  degrees  of 

freedom.  Used  for  asymmetric  loading. 

3  Isoparametric  quadratic  element  with  three 
nodes  per  side  and  R  and  Z  degrees  of  freedom. 

4  Element  to  be  used  on  mode  0  of  an 

asymmetric  loading  run  when  theta  loads  are 
used.  Torsion  in  connection  with  radial  and 
axial  loads. 

5  Element  to  be  used  on  mode  0  of  an 

asymmetric  loading  run  when  torsion  is  the  only 
load.  This  element  is  much  faster  then  type  4. 

9  Null  or  void  element.  This  type  must  be 
assigned  to  the  last  node  In  each  row  and 
column.  It  may  also  be  used  to  define  holes  in 
the  material.  It  is  automatically  assigned  to 
the  nodes  on  I  MAX  and  JMAX. 

L6  Material  number 


L7 

II  -if  input,  1 
increments  of  II 

will 

CO 

from 

11 

to  12 

I  n 

L8 

Jl  -  if  Input,  J 
Increments  of  Jl 

wl  1 1 

go 

from 

Jl 

to  J2 

in 
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Boundary  Condi t ions  -  The  flag  is  BOUN 


Each 
nodes . 

record  assigns  boundary  conditions 

to 

a  set  of 

LI 

11 

L2 

J1 

L3 

12 

L4 

J2 

L5 

R  boundary  condition  code* 

L6 

R  boundary  condition  value 
coef f i c lent 

or 

harmonic 

L7 

0  boundary  condition  code* 

LS 

0  boundary  condition  value 
coefficient 

or 

harmonic 

L9 

Z  boundary  condition  code* 

L10 

Z  boundary  condition  value 
coefficient 

or 

harmonic 

Lll-12 

not  used  at  present 

L13 

Sliding  boundary  condition  < 
only. 

;ode 

,  0  or  1 

L14 

Angle  in  degrees  of  the  line  the  node  must 
slide  on.  See  Figure  2. 

L15 

II  if  input,  1  will  go  from 
Increments  of  II 

11 

to  12  in 

L1C 

Ji  If  input,  J  wl 1 1  go  from 
increments  of  JI 

JI 

to  J2  in 

♦The  val id  codes  are: 

0  -  no  boundary  condition 

1  -  displacement  boundary  condition 

2  -  A  nodal  force  Is  to  be  applied.  Nodal  forces  are  the 

forces  on  the  entire  circumference. 
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R 

* 


Figure  2.  Sliding  Boundary  Condition 


Pressures  -  The  flag  is  PRES 

Each  record  puts  pressure  along  a  line  of  nodes.  P 
is  normal  to  the  element  surface  and  S  is  parallel  to 
it.  The  values  PR#  PTH  and  PZ  are  traction  loads  In 
the  coordinate  system  and  will  decompose  into 
pressure  and  shear  components  internally. 


LI 

11 

L2 

J1 

L3 

12 

L4 

J2 

L5 

Di rect i on 

0  if  the  pressure  is  to  be  applied  to  the 
elements  to  the  left  of  the  line  moving 
from  point  II,  J1  to  point  12,  J2. 

1  if  the  pressure  Is  to  be  applied  to  the 
elements  on  the  right. 

See  Figures  3  and  4. 

Note;  This  code  Is  correct  for  a 
right-handed  indexing  system.  For  a 
left-handed  indexing  system,  the  codes 
must  be  interchanged.  That  is,  0  is 
pressure  to  the  right  and  1  to  the  left. 
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L6 

PI 

Pressure  at  node 

11,  J1 

L7 

P2 

Pressure  at  node 

12,  J2 

L8 

SI 

Shear  at  node  11, 

J1 

L9 

S2 

Shear  at  node  12, 

J2 

L10 

PR1 

Radial  traction  at  node 

11, 

J1 

Lll 

PR2 

Radial  traction  at  node 

12, 

J2 

L12 

PTH1 

Theta  traction  at 

node 

11, 

J1 

L13 

PTH2 

Theta  traction  at 

node 

12, 

J2 

L14 

PZ1 

Axial  traction  at 

node 

11, 

J1 

LI  5 

PZ2 

Axial  traction  at 

node 

12, 

J2 
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X-N. 


Body  Force 

loads. 

-  The  Flag  Is 

BODY 

LI 

11 

L2 

J1 

L3 

12 

L4 

J2 

L5 

Radial 

(pwa ) 

body  force 

harmonic 

coef  f  1  dent 

L6 

Theta 

body  force  harmonic  coefficient  (F) 

L7 

Axial 

(F/L8) 

body  force 

harmonic 

coefficient 

Temperature 

l  Distribution  -  The 

Flag  Is  TEMP 

There  are 

two 

ways  in 

wh i ch  a 

temperature 

distribution  may  be  applied  in  this  program.  Only 
one  form  may  be  used  In  any  one  run. 


A. 


3. 


LI 

L2 

L3 


Constant  temperature  over  the  whole  body. 

The  first  record  must  be  the  word  CONSTANT 
enclosed  in  single  quotation  marks.  This 
record  must  be  followed  by  one  record  of  the 
following  form: 


LI  The  constant  temperature 
Radius  vs  temperature  tables 


The  first  record  must  be  the  character  string 
R  vs  T  enclosed  in  single  quotation  marks. 
This  record  will  be  followed  by  as  many  records 
of  the  following  form  as  necessary  to 
completely  define  the  temperature  of  each 
element : 


II  These  are  node  point 

J1  indices  which  define 


12 


the  area  over  which 
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L4 

J 

this  table 

appl ies. 

L5 

R1 

L6 

Tl 

Temperature 

at 

radius 

R1 

L7 

R2 

L8 

T2 

Temperature 

at 

radi us 

R2 

Material  Properties  -  The  flag  is  MATE 


This  program  accepts  isotropic  materials  only. 

For  nonstrain  or  temperature-dependent  materials# 
only  one  record  is  entered  for  each  material. 

For  materials  that  are  strain  or  temperature 
dependent#  several  records  will  be  entered  for  each 
material.  One  record  must  be  entered  for  each  entry 
in  the  table.  The  material  number  and  name  need  not 
be  entered  on  each  record  of  the  table.  If  the  table 
is  strain  dependent#  the  strain  value  in  the  first 
record  must  be  nonzero  and  the  table  must  be  so  input 
that  a  zero  strain  value  can  be  found.  If  the  table 
is  temperature  dependent#  the  temperature  value  in 
the  first  record  must  be  nonzero. 


LI 

Material  number 

L2 

Young's  modulus 
Kl* 

(E)  or  large 

deformat i on 

L3 

Poisson's  ratio 
K2* 

(v)  or  large 

deformat i on 

L4 

Alpha  for  thermal 

loading 

L5 

Temperature  for 
properties  apply 
dependent . 

which  these  material 
if  they  are  temperature 

L6 

Strain  at  which 

these  material 

properties 

apply  if  they  are  strain  dependent  (large 
deformation) . 


L7-9  A  24-character  material  name  enclosed  In 

quotation  marks. 


*See  Vol  II  for  definition. 
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NOTE:  Due  to  special  Internal  handling  of  materi.il 
properties#  da  not  use  L-Numbers  in  the  input. 
Always  specify  missing  parameters  by  commas. 

Selective  Print  Option  -  The  Flag  Is  SELE 

This  section  may  be  used  to  ask  for  output  of 
displacements#  stresses  and  material  properties  if 
the  overall  output  of  those  items  has  been  suppressed 
in  the  Output  section.  As  many  records  os  desired 
may  be  entered. 


LI 

11 

L2 

J1 

L3 

12 

L4 

J2 

Incremental  Loading  -  The  Flag  Is  INCR 

This  section  will  cause  the  program  to  iterate# 
applying  a  portion  of  the  load  each  time  and 
accumulating  the  results. 

LI  The  number  of  Increments  (N)  over  which 

the  loads  will  be  applied.  If  this  Is 
positive#  then  the  second  record  for  this 
set  must  be  entered  with  N  fractions  on 
it.  If  negative#  the  loads  will  be 
applied  in  an  arithmetic  progression. 

L2  If  a  nonzero  value  is  Input#  printout  of 

material  properties#  displacements  and 
stress/strains  will  be  made  each  Iteration 
for  those  nodes  specified  by  the  selective 
print.  The  complete  output  will  be 

printed  on  the  final  Increment  unless  the 
print  suppress  flag#  output  option  1#  Is 
turned  on. 

The  following  set  of  flags  tells  which  loads  will  be 
applied  incrementally.  If  the  flag  is  zero#  the  load 
will  be  applied  in  increments;  if  nonzero  the  full 
load#  if  any#  will  be  applied  In  each  iteration. 

L3  Pressure  and  traction  loads  PR#  PTH  and  PZ 
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L4 

Shear 

L5 

Temperature 

L6 

R  body  force 

L7 

Z  body  force 

L8 

R  nodal  load 

L9 

Z  nodal  load 

L10 

Nodal  displacements 

1 ncremental 

loading  fraction  record 

This  record  must  be  entered  If  N  Is  positive  and  may 
be  entered  to  override  calculated  values  if  N  Is 
negative. 

LI 

Fraction  for  1st  Increment 

L2 

Fraction  for  2nd  increment 

etc. 

for  N  (max.  of  100)  Increments 

iALSfi  Deformation  -  The  Flag  Is  LARG 

This  section  controls  the  large  deformation 
Iteration.  One  record  of  the  following  form  must  be 
entered. 

LI 

The  number  of  the  last  Iteration  to  be 
al lowed. 

L2 

If  nonzero,  printout  of  displacements  will 
be  made  each  Iteration  for  those  nodes 
specified  by  the  Selective  print  option. 

L3 

Convergence  factor.  Convergence  is 
attained  when  the  absolute  value  of  the 
largest  change  in  displacements  is  less 
chan  this  value. 

L4-103 

The  underrelaxation  factors  for  iterations 
1-99.  The  default  value  Is  1.0. 
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Geometry  Plots  -  The  Flag  is  GEOM 

This  section  is  used  to  request  the  plotting  of  the 
geometry,  either  deformed  or  undeformed.  Plots  for 
which  the  displacement  multiplier  is  zero  will  be 
produced  during  the  input  phase  at  the  time  the  plot 
section  is  encountered.  It  is  possible  to  plot  the 
geometry  at  any  point  in  its  development  by  putting 
geometry  plot  sections  in  at  the  desired  points,  such 
as  after  the  node  records,  etc. 


Note  that  data  will  not  carry  over  from  one  record  to 
the  next  in  this  section. 

LI  Axis  parameter  option 

1  if  the  R-axis  is  to  be  horizontal  to 
the  right 

2  if  the  Z-axis  is  to  be  horizontal  to 
the  right 


L2 


II 


These  indices  define  the  region  of 


L3 


J1 


the  geometry  to  be  plotted. 


L4  12 

L5  J2  J 

L6  The  displacement  multiplier.  If  zero,  the 

original  geometry  will  be  plotted.  If 
nonzero,  the  displacements  will  be 
multiplied  by  this  before  being  added  to 
the  node  coordinates  for  plotting. 


L7  X-axis  scale  factor* 


| 

/ 


♦The  X  and  Y  axes  here  refer  to  the  paper  rather  than  the 
geometry.  Not  all  of  these  parameters  need  be  entered.  The 
maximum  and  minimum  values  will  be  obtained  from  the  data  if 
not  specified.  Also,  both  the  scale  factors  and  the  paper 
size  need  not  be  specified.  The  relation  Scale  ■  paper  size 
/(Max  -  Min)  will  be  used  to  find  any  missing  data.  The  scale 
factor  will  always  be  used  if  entered.  If  the  scale  factors 
are  not  input,  they  will  be  calculated  to  keep  the  plot  on 
the  paper.  If  the  X  paper  size  and  scale  factor  are  not 
entered,  the  Y  scale  factor  will  be  used  for  both  axes. 
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L8 


Y-axis  scale  factor* 


L9  Minimum  value  on  the  X-axls* 

L10  Maximum  value  on  the  X-axis* 

Lll  Minimum  value  on  the  Y-axis* 

L12  Maximum  value  on  the  Y-axis* 

L13  X-axis  paper  size* 

L14  Y-axis  paper  size* 

L15  The  X  location  in  inches  of  the  beginning 

of  the  title. 

L16  The  Y  location  in  inches  of  the  beginning 

of  the  title. 

L17  The  letter  size  for  the  title  in  inches. 

L18-19  Not  used  at  present 

L20-25  Up  to  48  characters  of  plot  title 

information  as  an  alpha  string  enclosed  in 
quotation  marks. 

Accumulation  Control  -  The  Flag  is  accu 

This  section  must  be  entered  only  on  the  last  case  of 
an  asymmetric  loading  run.  Under  control  of  these 
records  the  displacements  and  strains  will  be  read  in 
for  each  mode  and  accumulated.  The  accumulated 
output  will  be  printed.  It  may  also  be  plotted  by 
putting  a  geometry  plot  request  after  the  desired 
displacements  have  been  calculated;  thus  there  may  be 
several  Accumulation  and  Plot  sections  alternating  in 
the  input. 


*The  X  and  Y  axes  here  refer  to  the  paper  rather  than  the 
geometry.  Not  all  of  these  parameters  need  be  entered.  The 
maximum  and  minimum  values  will  be  obtained  from  the  data  if 
not  specified.  Also,  both  the  scale  factors  and  the  paper 
size  need  not  be  specified.  The  relation  Scale  ■  paper  size 
/(Max  -  Min)  will  be  used  to  find  any  missing  data.  The  scale 
factor  will  always  be  used  if  entered,  if  the  scale  factors 
are  not  input,  they  will  be  calculated  to  keep  the  plot  on 
the  paper.  If  the  X  paper  size  and  scale  factor  are  not 
entered,  the  Y  scale  factor  will  be  used  for  both  axes. 


Note  that  to  plot  a  transverse  cut#  I  goes  from  1  to 
the  maximum  number  of  nodes  In  the  Radial  direction 
and  J  goes  from  1  to  19. 

LI  Option  Flag 

1  Output  will  be  an  axial  cut  at  angled. 

2  Output  will  be  a  transverse  cut  at  the 
axial  di stance  Z. 

L2  e  or  Z  depending  on  LI. 

Stabi 1 i tv  -  The  flag  is  STAB 

This  section  Is  used  to  input  the  data  for  the 
stability  calculation.  Run  option  20  must  be 
specified  for  the  calculations  to  be  made.  The  stack 
is  assumed  to  be  a  vertical  stack  beginning  and 
ending  with  a  rubber  pad;  hence  the  number  of  shims 
is  assumed  to  be  one  less  than  the  number  of  rubber 
pads. 

LI  The  thickness  of  each  rubber  pad. 

L2  The  thickness  of  each  shim. 

L3  The  inside  diameter  of  the  stack  if  It  is 

cyl Indr ical . 

L4  The  outside  diameter. 

L5  The  number  of  rubber  pads. 

L6  The  Young's  modulus  for  the  rubber  pads. 

L7  The  convergence  accuracy  factor.  When  (1 

-  old  buckling  load/new  buckling  load)  is 
less  than  this  value#  the  iteration  has 
converged. 

L8  The  maximum  number  of  iterations  to  be 

al lowed. 

End  of  Harmon i c  -  The  flag  is  END 

An  end  record  must  terminate  the  input  for  each 
portion  of  the  run.  This  signals  the  program  to  stop 
reading  and  solve  the  problem  as  defined. 


44 


SAMPLE  1MEULI 


Two  cases  are  shown  here.  They  are  illustrations  only  and  do 
not  necessarily  reflect  real  problems. 

Sample  case  1  is  a  flat  bearing.  It  is  set  up  in  a 
right-handed  system.  Figure  5  Is  a  plot  of  the  geometry,  and 
arrows  have  been  drawn  in  to  show  the  R,  Z,  I  and  J  axes. 
This  Is  purely  a  bearing  with  no  additional  parts  added.  The 
output  for  this  case  is  not  given  in  this  document  since  its 
form  is  similar  to  that  of  case  2. 

Sample  case  2  is  a  spherical  bearing  with  a  steel  attachment 
on  one  end.  It  is  oriented  in  a  left-handed  system.  Figure  6 
is  a  plot  of  the  geometry,  and  arrows  have  been  drawn  in  to 
show  the  axes.  Note  how  the  input  sections  are  ordered  to 
achieve  the  desired  results.  The  output  from  this  case  is 
given  here  to  show  Its  form. 
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SAMPLE  PROBLEM  1 


SIMPLE  PROBLEM  1  ; 

GENERAL  DATA  2 

l#  i»  6»  14  ; 

RUN  ; 

3#  4  ; 

OUTPUT  ; 

6#  9#  14  : 

BEARING  2 

l*  6#  l#  18#  1*  1* »  C >  C #  3.C#  3 .  G  »  0*  2.0#  0#  2.C# 

.55#  2#  2#  .3#  3#  3  2 
BCUNOARY  CONO  2 

1#  1#  1#  14#  1#  0  2 
1#  1#  6#  1#  0#  G#  0#  C#  1#  0  2 
BOCY  FORCE  2 

1#  l#  6#  4#  0#  C»  1C.  2 
1#  6#  6#  9  2 
1#  11#  6#  14  2 
1#  4#  6#  fc#  G#  0#  3C.  2 
1#  9#  6#  11  2 
MATERIALS  2 

1#  1.0E6#  .3###*  •  STEEL  SHIM*  2 
2#  3C0.#  .^SSS9,#,,  *RU8Pf R  PAC*  2 
GEOMETRY  PICT  2 

1#  l#  l»  fc*  14#  0#  3.*  3.#  -.333»»  -.333####  l.C#  3.5#  .12# 
L20  'SAMPLE  1*  2 

1#  l»  1#  t#  14#  5.0#  4.#  4.#  C##  C#»»#  1.0#  P.5#  .12# 

L20#  'OE  FORMEO  SAMPLE  1*  2 
ENC  2 


SAMPLE  PROBLEM  2 


SAMPLE  PR06LEP  2*  ZERCTH  HAkHCNIC  i 
GENERAL  CAT  A  i 
I#  l#  6#  IE#  111#  0#  1  • 

RUN  ; 

3#  4#  9. 

NCCES  f 

I#  1#  0*0#  9*0  • 

6#  1#  3.0#  9.0  ; 

BEARING  ; 

1#  6#  3#  LB#  2#  C#  IC.C#  C#  C#  2.C#  3.C#  10. C#  12. C#  9.0#  II. 0# 
•3#  2#  2#  .24#  3#  3* 

LINES  ; 
l#  l#  b#  l  • 

l#  l#  l.  3  i 

6#  1#  b#  3  # 

GRIO  S 

I »  I  *  6#  3  J 
TYPE  ; 

l  #  1 #  5#  2 »  2#  l# 

BOUNDARIES  S 
1#  2#  1#  It#  I#  0  J 
I#  l#  6#  1#  1#  0#  l»C#  1#  C  » 

PRESSURE  S 

I#  16#  b#  lb#  1#  50.#  !0.  » 

RATER  I ALS  i 

I#  1.CE6#  .3### #  *  STEEL  '  ; 

2#  300.#  .49999####  *  RL3  6£R  *  i 
GEOMETRY  PLOT  S 

2#  l#  l#  b*  lb#  C*  2.0*  2.C#  B.C»*  -l.C#*##  2.C#  d.3*  .12  * 

L2C#  •  SAPPLE  2  LN Cfc FCRMtC  *  5 
ENC  OF  HARMONIC  i 

SAMPLE  PRC6LEM  2#  FIRST  MRPCNIC  ? 

GENERAL  DATA  i 
I#  I#  6#  lb#  111#  1#  1  • 

RUN  ; 

3#  4#  9# 

NCCES  * 

1#  1#  C.O#  9.C  i 
6#  1#  3.0#  9.C  • 

BEARING  ; 

1#  b#  3#  L 8#  2#  0#  IC.C#  C#  C»  3.C#  3.C#  10. C#  12. C#  9.0#  11.0# 
•  3#  2#  2#  .24#  3#  3# 

LINES  i 
1#  1#  b#  I  S 
I#  i*  l#  3  l 
b#  i#  b#  3  # 
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€1110  ; 

1*  If  6*  3  5 

Tm  ; 

1#  I#  3*  2f  2#  l  5 
BOUNDARIES  5 

It  1*  6 1  1*  1#  0 1  l»0i  l#  0  5 
PRESSURE  5 

tt  1«  6f  16»  C*  100.  f  100.5 
MATERIALS  5 

It  1.0E6f  .3 tttt  •  STEEL  •  5 
it  300. »  .49 <Mtttt  •RUBBER'  5 
ACCUMULATION  5 
1«  0  I 

GEOMETRY  PLOT  5 

it  If  If  6f  16#  Sf  2 .Of  2.Cf  8. Cf#  -1,0.#*.  2. C»  8,3#  *12  * 
L20#  •  SAAPIE  2  OEFORPED  S»*S 
ENC  OF  HARMONIC  5 
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SAMPLE  OUTPUT 


The  output  shown  is  for  Sample  Input  Case  2.  Some  of  the 
pages  are  combined  to  save  space,  and  in  some  cases  the  full 
output  is  not  shown.  Note  that  the  input  for  the  second  pass 
has  been  heavily  abbreviated  since  most  of  it  is  the  same  as 
for  the  first  pass. 

Figure  7  is  a  plot  of  the  deformed  geometry  with  the 
displacements  multiplied  by  5. 


SAMPLE  2  DEFORMED  5X 


Figure  7.  Deformed  Sample  2 
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APPENDIX  A 


Ssrvlse  Life  Program 


This  program  is  designated  as  S3359SL.  It  is  designed  to 
accept  a  tape  or  tapes  created  by  S3359,  when  output  option  20 
has  been  specif ied,  and  accumulate  a  selected  set  of  stresses 
and  strains. 

The  energy  is  calculated  for  each  element  from  the  principal 
stresses  and  strains  as  the  accumulation  is  done. 

The  accumulated  energy,  stresses  and  strains  and  the  principal 
stresses  and  strains  calculated  from  the  accumulated  stresses 
and  strains  are  printed  out  at  the  end  of  the  run.  The 
locations  of  the  highest  energy  and  the  highest  value  of  the 
principal  stresses  are  printed  out  also. 

lUCLUl 

When  the  S3359  runs  were  made,  a  pair  of  integers  giving  the 
date  and  time  for  each  set  of  data  put  on  the  tape  was  written 
out.  To  select  the  desired  sets  of  data  from  the  tape  for 
accumulation,  the  date  and  time  of  each  of  the  desired  sets 
must  be  input  exactl v  as  they  were  written  out  and  In  the  same 
order  as  they  appear  on  the  tape. 

The  input  is  via  FREFRM  as  follows: 


Title  -  A  title  record  that  will  be  printed  on  the  top  of  each 
page  of  the  output.  This  must  be  only  one  card. 

Control  Records 

As  many  records  as  desired  to  select  the  data  sets. 

LI  The  data  as  printed  by  S3359. 

L2  The  time  as  printed  by  S3359. 

L3  The  factor  by  which  the  stresses  and 

strains  will  be  multiplied.  The  default 
value  is  1.0. 

An  end  of  file  stops  the  reading  and  accumulating  and  causes 
the  final  output  to  be  produced. 
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APPENDIX  B 


S3359F  Fourier  Coefficient  Generator 


This  program  computes  the  Fourier  series  coefficients  by  use 
of  Trapezoidal  Integration  Formula#  then  fits  the  function 
f  ( 0 )  with  a  Fourier  series  curve  fit.  The  output  Includes  a 
set  of  cards  containing  the  Fourier  series  coefficients  In  a 
form  compatible  with  the  Input  requirements  of  program  S3359. 

The  function  to  be  fit  with  a  Fourier  series  curve  fit  must  be 
either  odd  or  even  having  period  2n.  The  Integration  for  the 
coefficients  is  performed  over  the  Interva  1  [o,«]  and  then 
multiplied  by  two. 

Preparation  oL  input 

The  input  is  read  In  in  free  form.  (The  rules  of  free-form 
Input  are  given  in  the  S3359  document  in  detail.)  The  first 
card  is  a  title  card.  The  second  card  contains  data 
specifying  the  symmetry  of  the  function#  the  error  tolerance 
which  is  compared  with  the  least-squares  error  approximation# 
the  maximum  harmonic  which  is  to  be  calculated#  and  a  plot 
flag  specifying  what  is  to  be  plotted.  The  next  set  of  cards 
contains  values  of  e  and  f<0)  where  0  Is  in  degrees. 


Title  Record 

May  contain  any  desired  alphanumeric  information 
Control  Record 

I OEF  IOEF-0  if  f ( 0 )  Is  even,  1  If  f(0)  is  odd 

I SYM  Specifies  the  symmetry  of  f(0)  as  shown  In 

Tables  B-l  and  B-2 

E  Real  number  error  tolerance  allowed  for 

Fourier  series  curve  fit 

N  Integer  number  specifying  the  maximum 

harmonic  to  be  computed.  Maximum  of 
fifteen. 

I  PLOT  I PLOT-O#  no  plot 


73 


IPLOT-1,  plot  all  harmonics 

IPLOT-2,  plot  final  harmonic  only 

IPLOT-3,  plot  specified  harmonics 

IH(15)  Integer  number  specifying  harmonic 

numbers  to  be  used  in  plotting. 
Maximum  of  fifteen.  Required  only 
for  I  PLOT-3. 

Input  Data  Cards 

0  Real  value  of  6  in  degrees 

f(e)  Real  value  of  f(e) 

The  title  card  must  be  ended  with  a  semicolon.  All  other  data 
is  separated  by  commas  and  ended  with  a  semicolon.  In 
choosing  N  and  IH,  the  user  must  keep  in  mind  what  type  of 
symmetry  (IOEF  and  ISYM)  he  Is  working  v/ith.  For  example,  if 
I0EF  ■  0,  and  ISYM  ■  1,  only  even  harmonics  are  calculated 

since  the  odd  harmonics  are  zero.  Hence,  N  and  the  values  of 
IH  should  be  even. 

If  the  Fourier  coefficients  are  to  be  computed  for  more  than 
one  function,  the  data  for  each  function  Is  preceded  by  a 
title  record. 

Many  functions  of  interest  contain  some  type  of  symmetry.  In 
these  cases  certain  terms  in  the  Fourier  series  will  be  zero 
and  therefore  are  not  calculated.  The  function  f,  in  this 
case,  is  always  either  odd  or  even.  The  following  rules  of 
symmetry  apply. 
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EVEN  FUNCTION 
Table  B-l 


Code 

Function 

Fourier  Terms 

ISYM=0 

f(  e)-f(-e) 

no  sine  terms 

ISYM-1  or  2 

f(0)-f(l8O+e)-f(i8O-e) 

no  sine  terns; 
no  odd  cosine 
terms 

1 SYM=3 

f(0)-f(l8O+e)--f(i8O-e) 

000  FUNCTION 

Table  B-2 

no  constant 
term 

no  even  cosine 
terms 

no  sine  terms 

Code 

Function 

F.fl.ur  Lar.  Iaroa. 

1 SYM-0 

f(-e)  — f  (e> 

no  constant 
term 

no  cosine  terms 

ISYM-1 

f(e)-f(i80+e)«-f(i80-e) 

no  constant 
terms 

no  cosine  terms 

no  odd  sine 

terms 

ISYM-2or3 

f ( ©) -f C 18  0- 0) »- f(180+e) 

no  constant 
term 

no  cosine  terms 

no  even  sine 

terms 

Note,  for  example,  in  the 
ISYM-1,  that  though  only 
calculated,  some  of  these 
most  likely  will  not  be 


case  where  the  function  is  odd  and 
the  even  sin  terms  are  being 
might  also  be  zero.  These  numbers 
zero  in  the  output  but  only  small 
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numbers.  The  user  must  decide  whether  or  not  they  are  zero 


Output 

The  entire  input  data  is  listed  giving  the  FORTRAN  name.  The 
computed  Fourier  coefficients  are  written  with  a  least-squares 
error  for  the  curve  fit  up  to  that  particular  harmonic.  A 
deck  of  cards  containing  the  Fourier  series  coefficients  is 
punched.  The  cards  are  punched  in  a  form  compatible  with  the 
input  requirements  of  Program  S3359. 

If  a  plot  Is  specified,  a  plot  of  the  Fourier  series  curve  fit 
is  given  on  12-Inch  grid  paper.  The  actual  curve  is  marked  by 
X's,  and  the  Fourier  series  curve  fit  is  drawn  with  a  solid 
line. 


Sample  Input 
PRESSURE  LOAD  1; 

1.1. . 001.12.3.2.12; 
0  • ,  2 . ; 

45. . 2.; 

90. . 2.; 

90,-2; 

135.,  -2; 

180.,  -2; 
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APPENDIX  C 


370  OPERATING  INSTRUCTIONS 


PROGRAM:  STRESS  ANALYSIS  OF  AN  AXISYMMETRIC  BODY 

WITH  ASYMMETRIC  LOADS 


TAPES 


)ens  i  ty 

2 

5 

a 

10 

L 2 

s 

a 

ie 

2  ! 

j  a 

ie 

2 

5 

8 

lfl 

2 

5 

a 

is 

)ata  Set 

Name 

- 

SERVLIFE 

1 nput/ 

Output 

I 

(i) 

xl 

0 

I 

1  o 

I 

!  0 

xl 

0 

Reserve 

List 

Scratcl1 

L  1  SC 

R  II  be 

lIsc 

l|sc 

R_ 

lx-  BC 

- 

CARD  OUT 

YES 

1 

(no) 

PLOTTER 

YES 

1 

W 

ERROR  PROCEDURE: 

Flush 

ESTIMATED  MAXIMUM  RUNNING  TIME!  See  user  request 
SPECIAL  INSTRUCTIONS:  None 
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